The vibrational density of states of a disordered gel model 
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I mixture of colloidal patchy particles with two and three patches for different relative compositions 

^ X2- Ai low temperature, this system forms a thermo-reversible gel, i.e. a fully bonded network 

' of chains of two-patches particles, in which the branching points are provided by three-patches 

CN particles. For all the compositions we find in the vDOS a pronounced peak at low frequency whose 
> 

^2 height grows on increasing the fraction of two-functional particles, or equivalently with the average 

^ length of the chains. To identify the various spectral features, we compare the vDOS of the whole 

^ system with the one of small representative structures of the network and with the vDOS of a long 

^ linear chain of two-patches particles and we find that these structures are indeed able to rationalize 

>• the various peaks in the vDOS of the full system. At large X2 the vDOS of the gel and of the long 

^ chain show remarkable similarities. Analyzing the dispersion relations and the spectrum of the 

5-! 



linear chain we show that the excess of low frequency modes, the analog of the boson peak in 
glassy disordered systems, arises from the strong coupling between rotations and translations. 
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I. INTRODUCTION 



The study of the vibrational density of states (vDOS) is central in the investigation of 
solids [HE], since, apart from the structure, it is probably the most relevant and direct quan- 
tity to characterize a system. E.g., within the harmonic approximation, the vDOS (which 
we will denote by g{uj), where u is the frequency), allows one to calculate the temperature 
dependence of the specific heat, thermal conductivity, as well as the free energy. It has been 
known for a long time that the vDOS of disordered systems shows anomalies, i.e. has fea- 
tures which are not observed in crystalline systems. A well-known example is the vDOS of 
fractal aggregates for which one finds that g{u) shows at low u a power-law dependence with 
a fractal exponent, instead of the typical Debye law w^^^, where d is the dimensionality of 
space [3H3] . Deviations from the Debye law are also often observed in atomic and molecular 
glasses and have been in the focus of interest of many studies [5H2U]. The most remarkable 
anomaly is an excess in the number of modes, as compared to the Debye prediction, at low 
frequencies (i.e. for an atomic system in the THz range), a feature which is commonly named 
"boson peak" and which is responsible for the anomalous behavior in the specific heat or 
thermal conductivity of glasses [9l [121 [131 [IS] . Although some evidence has been provided for 
the existence of a correlation between the fragility of a glass-forming system and the strength 
of its boson peak, other results suggest that no such correlation exists [3 |2I1 [22] • It has 
been proposed that this correlation has a structural origin [20l |21] , since strong glass-formers 
have an open network structure (due to the presence of directional/ionic bonds) and tend to 
have a strong boson peak, while fragile liquids have a structure that is similar to the one of 
the hard sphere system, i.e. van der Waals like interactions, and usually a weak boson peak. 
Similar results have been found in the Sea;Asi_a; system for which neutron scattering shows 
that the number of low frequency modes is related to the average coordination number of 
the atoms |23j . 

While much effort has been devoted to the investigation of dense glasses, much less is 
know for the case of thermoreversible gels, i.e. disordered arrested states at low densities |21]. 
One of the reasons for this is that the realistic modeling of thermoreversible gels has become 
possible only recently [25143T] . Models include patchy particles of low valence, particles 
interacting via competing short and long range interactions, three-body potentials, dipolar 
particles, and others [21]. In all these models, the interaction potential makes that the 
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particles have only a small coordination number which in turn allows to generate very open 
disordered structures, composed by chains of bi-coordinated particles joined by a small 
number of three- (or more) coordinated particles which act as junctions in the network. 
In the last years the structural and relaxational properties of these thermoreversible gel 
models have been extensively investigated in the attempt to deepen our understanding of the 
mechanisms responsible for the slowing down of the dynamics with decreasing temperature 
and to clarify the differences between gels and glasses |2^ . 

In the present work we extend this line of investigations and focus on the vDOS of one 
of these models, a binary mixture of particles with, respectively, two and three attractive 
patches. By changing the relative composition of the mixture we can change in a controlled 
way the structural properties of the system ^32j and, in particular, the length of the bi- 
functional chains connecting the three-functional particles. This flexibility will allow us to 
study how a change of the structure influences the properties of the vDOS as well as the 
behavior of the system at low frequencies, where a strong coupling between rotational and 
translational degrees of freedom is found. 

II. MODEL AND COMPUTATIONAL DETAILS 

The model used in the present study has been presented in detail in Ref. |32j. Hence 
we provide here only the necessary amount of information required to make the article 
self-contained. 

A. The interaction potential 

The particles are modeled as spherical rigid bodies of diameter a and mass m. The 
surface of particle i is decorated with Mj sticky spots called patches. These patches are 
located at a distance a/2 from the particle's center of mass. The total interaction between 
particles i and j is 

V{t,j) = VcKS^j) + Vp{t,3), (1) 

where Vcm is the interaction between the centers of mass of the two particles, while Vp 
is the potential acting between the patches. So the total potential of the system is Vtot = 
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and 
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Vp{i,j) = -e^^exp 



(3) 



p=l q=l 



where Vij is the distance between the centers of mass of particles i and j and r^J is the 
distance between patch p on particle i and patch q on particle j. With this choice of 
parameters, the resulting potential resembles a hard sphere potential complemented by an 
attractive square well-like potential if patches on different particles are close to each other. 
In addition, the very short range of Vp guarantees that the so-called "single bond per patch 
condition" is fulfilled. This condition means that the maximum number of possible bonds 
that particle i can form is equal to Mj, and it allows to make an analytical description of 
the structural and dynamical properties of the model [271 1331 [33]- By choosing e = 1.001, 
one finds that the depth of the V{i,j) potential is Uq = 1.0. 

B. Details of the simulations 

The studied system is a binary mixtures of N2 bi-functional particles, Mj = 2, and 
N3 three functional particles, Mj = 3, with = A^2 + ^3 = 1000 at packing fraction 
(j) = l^y^ = 0.1, where V is the total volume of the system. Patches are arranged on 
the poles for particles with two patches (bi-functional particles) and equally spaced on the 
equator for particles with three patches (three functional particles). We have studied six 
different systems having the following concentrations of bi-functional particles X2 = N2/N: 
0.2, 0.5, 0.66, 0.75, 0.85, and 0.9. This gives us the possibility to investigate how the average 
valence M = 3 — X2 influences the structure and the vibrational properties of the system. In 
the following, energy is measured in units of Mq, distance in units of a, and time in units of 
yjmo'^ je. Temperature, T, is measured in units of energy, setting the Boltzmann constant 

The configurations we analyze were generated via Brownian dynamics [32] by equilibrat- 
ing the systems at T = 0.055 (a T at which essentially most of the possible bonds have 
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already been formed) and by quenching the final equilibrium configurations to T = 0.04. 
The simulation at T = 0.04 is performed until the few remaining isolated monomers attach 
to the spanning cluster. Indeed, at such a low temperature, the probabihty of breaking a 
bond, proportional to e"^'''^^^, is of order 10^^^ and hence no bond breaking occurs during 
the simulation (~ 10*^ time steps, where each time step correspond to 5t = 0.001). We 
stress that, in the resulting structure, all the particles belong to the same cluster and the 
connectivity of the system is not altered by the dynamics, i.e. in the following we study 
only the vibrational dynamics and not relaxation and thus the system can be considered as 
a chemical gel. Note that this vibrational dynamics is not harmonic since the system is able 
to overcome local barriers, but which are not associated to bond breaking events. 



C. Method to evaluate the harmonic density of states 

We have computed the normal modes and eigenfrequencies of the system by determining 
the properties of the local potential energy minimum configuration q*, i.e. a configuration for 
which VqVtot(q)|q, = 0, where q are the generalized coordinates (center of mass positions 
and orientational angles). This local minimum was obtained by taking a configuration of 
the system as a starting point of a minimization procedure in the potential energy Vtot 
and the so obtained local minimum is q*. (Note that this minimum configuration is often 
also called "inherent structure" (IS) [35l|36].) In practice we have used a conjugate gradient 
algorithm [37] which is here a fast procedure since the system needs only to relax the thermal 
fiuctuations without altering the connectivity of the system. 

Within the harmonic approximation, the local minimum of Vtot can be written as a bilinear 
form in qj and q^ 



i=i j=i 

where G is the number of degrees of freedom and 



Hi, = 



d'Vtotici) 



dqidqj 



(5) 



The eigenvectors of the matrix Hjj give the normal modes of our system and the square 
root of the eigenvalues give the eigenfrequencies. The normalized distribution of the latter 



is thus the vDOS g^u). 

To compute the Hessian Hjj at the IS we have used a finite difference method. Since 
bifunctional particles have patches located on the poles, only two angles are needed to iden- 
tify the direction of the axis that joins them. This lowers the total number of degrees of 
freedom from 6A^ — 3 to 5N2 + QN3 — 3. The subsequent diagonalization was performed using 
the LAPACK package. To improve the statistics we have averaged over 20 — 50 configura- 
tions. Although these configurations are not independent (since they share the same network 
topology), this averaging significantly smooths the vDOS. We associate such improvement 
to sampling of the several secondary minima composing the same metabasin |38j . 



III. RESULTS 



A. Structure of the investigated systems 

The networks studied in this article are composed by chains of bi-functional particles 
that are cross-linked by the three functional particles, acting as hinges between the chains. 
Since bonds do not break during the time scale of the simulation, the system behaves as 
an almost fully connected chemical gel [32]. The possibility of tuning the concentration of 
junctions allows us to generate structure that have different chain lengths. Specifically, the 
average chain length, i.e. the average number of bi-functional particles between two three 
functional ones (a length that is related to the average distance between branching points 
of the network) is given by [32] 

For the smallest value of X2 studied here we have thus I ~ 1.2, while for the largest we 
obtain / = 7.0. The ability to predict theoretically / is made possible by the fact that the 
structure is essentially generated from an equilibrium configuration and hence can be quite 
accurately described by the theory of Flory and Stockmayer[33l 139] . 

To provide a visual representation of the networks generated for different values of the 
mean connectivity M, we show in Fig. [ijtwo snapshots of the systems with X2 = 0.2 and 
X2 = 0.9, where the most mobile particles are highlighted |40j. From the figure we recognize 
that for X2 small, i.e. high connectivity, the system seems to be more heterogeneous than 



6 



FIG. 1. A snapshot of the system at T = 0.04 for X2 = 0.2 (a) and X2 = 0.9 (b). Small spheres: 
Bi-and three-functional particles with low mobility (red and blue, respectively). Large spheres: Bl- 
and three- functional particles with high mobility (light and dark green, respectively). A particle 
is defined to be mobile if it has a mean squared displacement > 0.5 for X2 = 0.2 and > 7.0 for 
X2 = 0.9. 

if X2 is large. This is due to the fact that in the former case the coexistence curve is closer 
to the investigated state point and thus fluctuations are larger [52] • Furthermore we see 
that for small average connectivity the most mobile particles form chains whereas for small 
X2 their spatial distribution seems to be rather random. It is evident that the chain-like 
structure found at low connectivity will give rise to low frequency modes in the vDOS, and 
below we will discuss this in more detail. 
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FIG. 2. Density of states 5(0;) for different fractions of bi- functional particles X2. At the lowest 
X2 there is one peak and a long tail. With increasing X2 a new peak appears while the minimum 
between the peaks becomes more pronounced but remains at the same position [yimin — 0.38). 
Consecutive curves are shifted vertically by 0.5. 

B. Density of states 

In this section we discuss the density of states, calculated within the harmonic approx- 
imation, for different concentration of bi-functional particles and correlate specific features 
of g{uj) with specific vibrational modes of the system. 

The harmonic density of states for several values of M is shown in Fig. [2J For small X2 
the vDOS has one main peak at small frequencies and a tail that extents to high u. With 
intermediate and large X2 the vDOS displays two peaks, one located for a; ^ 0.38 and one 
centered around w ~ 1, followed by a broad tail extending up to w ~ 10. Apart from the 
case X2 = 0.2 (corresponding to / ~ 1.2), which shows a clear decrease in the number of 
low frequency eigenmodes, no significant dependence on X2 is observed in the fraction of 
modes contributing to the two main peaks and to the tail of the vDOS. This suggests that 
increasing the length of the chains does not strongly influence the total number of modes in 
the vDOS at low frequencies, but only the shape of g{uj). Indeed, with increasing X2, the 
low frequency peak moves to lower frequencies and increases its height while the peak at 
intermediate frequencies becomes better resolved and a shoulder appears at u ~ 0.6. For 
completeness, we note that a shoulder is also present in the low frequency peak {u ^ 0.1), 
which turns into a minor peak when X2 = 0.9. 

Note that the peak at low u is observed at frequencies that are significantly smaller than 
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the typical frequencies of the model, which are on the order of unity (due to the choice of 
the units). Since the mentioned peak occurs at frequencies that are by a factor of 3-10 times 
smaller, we can tentatively identify this peak as a boson peak. From the figure we thus 
can conclude that there is a clear correlation between the structure of the system (open vs. 
compact network) and the height of the peak (high/low). 

To better understand the nature of the oscillations in the different regions of the spectrum, 
we study the participation ratio, Pr, of each mode i. This quantity is defined as 

Priuj)=(^N^Je-f^ , (7) 

where e" is the displacement vector of the n-th particle of the i-th normal mode. The 
participation ratio is a standard measure of localization of a mode[l2l H3]. For example, for 
a translation of the system Pr{u)) = 1, while for a vibration of a single particle Pr{u}) = 1/N. 

In dense systems the direction of the vibrations of a given particle, i.e. e", is basically 
independent of the direction of the bonds that the particle forms with its neighbors. For a 
system that has long chains, as the one studied here, this is no longer the case and hence 
it is of interest to be able to discriminate between longitudinal and transverse modes with 
respect to the bond direction (i.e. along the local curvilinear direction of the chain). To 
this aim we quantify to what extent the displacement of the i-th normal mode is parallel to 
the local direction of the bond. By writing the direction joining bonded particles n and j 
as v" = Vnj/rnj, the observable 0{u) is thus defined as 

1 ^ 

n=l 

where j is an index running over the bonded neighbors of particle n. Since v" as well as e" 
are normalized, 0{uj) G [0, 1]. In the case of a uni dimensional chain, a longitudinal mode 
has 0{u)) = 1 while a transverse one has 0{u)) = 0. 

The functions Pri(^) and 0{uj) are shown in Fig. [sjfor X2 = 0.2 and X2 = 0.9, together 
with the corresponding vDOS for an extended frequency range. On increasing X2 both 
functions become more structured, suggesting that different frequency regions are populated 
by modes with different characteristics. We start by discussing the X2 = 0.9 case. The data 
for 0{u)) shows that modes for u > 1.5, i.e. modes that are not in the two peaks seen in 
g{uj), have 0{u) ~ 1, suggesting that such modes are mostly vibrations along the bond 
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FIG. 3. g{u}) (top), Pr(w) (middle), and 0{uj) (bottom) for X2 = 0.2 and X2 = 0.9. 

directions (reminiscent of longitudinal modes). In contrast to this, modes with u < 1.5 are 
characterized by a rather small value of 0{uj), thus describing oscillations orthogonal to the 
chain direction, i.e. are similar to transverse modes. Note that the transition of 0{uj) ^ 1 to 
values < 0.5 at ~ 1.5 is quite sharp, indicating a well defined cross-over from longitudinal 
modes at high frequencies to transverse modes at low u. The participation ratio, which 
never exceeds the value 0.4, suggests that all modes are mostly localized. We recall that the 
value usually taken as indication of a collective oscillation, i.e. of propagating phonons|i3]. 
is Pr ^ 0.5. Note that such propagating phonons are not even seen at low frequencies where 
one does expect acoustic modes. The absence of such excitations indicates that on the length 
scale of the simulation box the system is still very heterogeneous and as a consequence the 
modes cannot propagate. However, it can be expected that for larger simulation boxes such 
propagating modes will indeed be present and the rise of Pr{(jj) at small co can be taken as 
indication that one is approaching this acoustic regime. Also at values of u in which g{u) 
shows the peak around ~ 1.0 we see that Pr{uj) is significantly larger than zero. Thus we 
conclude that also these modes show a significant delocalization. As we will show below, 
these modes correspond in fact to excitations in which a chain makes transverse oscillations, 
in agreement with the small values of 0{uj) in this w— range. Last but not least we see that 
Pr{uj) is also relatively large if u > 2.0. This is (probably) related to the fact that at elevated 
frequencies the modes are a hybridization of (relatively) localized modes, thus should not 
really be considered as extended objects |13]. 
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FIG. 4. Studied network elements: a) A curved chain formed by 18 particles; b) Three chains 
linked together, 18 particles in total; c) A closed ring. All these structures have been found in 
X2 = 0.9 configurations at T = 0.04. 

In the case of X2 = 0.2, 0{u) never reaches very small values and the difference between 
the region of the peaks and the tail is less pronounced. This indicates that even modes 
with u < 3 cannot be described to be transverse (with respect to the bond direction). 
This is consistent with the very short average chain length observed for small values of 
X2- Furthermore we see that Pr{u}) becomes relatively large as soon as w < 1.0, i.e. for 
these frequencies the modes are significantly more collective than for the case X2 = 0.9. 
This is in agreement with the fact that for X2 = 0.2 the structure of the system involves 
smaller length scales than the ones for large values of X2 and hence acoustic like modes can 
propagate even at relatively high frequencies. From 0{uj) we also recognize that at large 
frequencies the modes appear to describe motion parallel to the bonds, i.e. 0{uj) is large. 
Since the participation ratio is relatively small, this suggests that these longitudinal modes 
are localized. This feature is thus directly related to the structure of the system which, for 
small X2, shows short rings of bonded particles, as displayed in Fig. [l^. 

To gain insight into the nature of the modes in the different regions of the spectrum 
we compare the vDOS of the system with the vDOS of isolated elements of the network, 
found in the simulated configurations. In particular we focus on three different structures: 
A bended chain of bi-functional particles, a star configuration in which three chains of bi- 
functional particles are connected via a central three-functional particle and a ring of bl- 
and three-functional particles. The first two structures are typical elements for the network 
with large X2, while the ring is relevant of the case with small X2- Figure |4] shows the three 

11 



3.0- 



1.0- 





1 .... 1 ... . 

_ bulk = 0.9 




— chain 




■ ■ ■ star 




- ■ ring 


1 1 * 1" 1 1 1 1 





J I I I I I I I I I I I I I I I I I l_ 



* 0.5 -m^o ° *^*^^^ ®> * * 

rt°o » * 




FIG. 5. g{uj) (top), Pr{uj) (middle), and 0{uj) (bottom) for the structmes pictured in Fig. |4j a ring 
(filled diamonds), a star of chains (stars) and a chain (open circles). For reference we include also 
the g{uj), Pr{u}) and 0{uj) of the entire system for X2 = 0.9. 



selected structures, composed respectively of eighteen, eighteen and seven particles. The 
large number of particles in the first two structures allows us to explore a larger and better 
resolved range of frequencies. For the chain the number of degrees of freedom, and therefore 
the number of eigenvalues obtained, is 92; for the star structure it is 94, and for the ring 38. 
We have also computed the density of states for chains whose length is comparable to the 
average chain length (/ ~ 7) and we have obtained the same qualitative results [H]. 

For each of these structures we have computed the density of states after the usual 
minimization stage. The results are shown in Fig. |5j The chain and the star configurations 
have density of states that are remarkably similar to the one of the system with X2 = 0.9. 
Instead, the ring's g{oj), is qualitatively more similar to the density of states of the system 
with X2 = 0.2, i.e. it does not have a pronounced peak at a; ~ 1.0. Figure [5] also shows the 
participation ratio and the 0(w) for the selected structures. Some of the features found in 
the bulk system are also observed in the eigenmodes of the examined structures, in that for 
u > 1.5, 0{uj) ~ 1 for the chain and the star while 0{u) ~ 0.2 for the ring. Moreover, in 
the frequency region u < 0.25, both 0{uj) as well as Pr{u)) are larger for the ring than for the 
chain and the star. Again these features are very similar to the ones found when analyzing 
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a) Long wavelength transverse acoustic mode 




c) Longitudinal acoustic mode 




w = 2.0 




d) Longitudinal optical-lilce mode a-*^ 




FIG. 6 



Four characteristic vibrational modes of the bended chain. 



the X2 = 0.2 and X2 = 0.9 systems. Hence, despite the small number of degrees of freedom of 
these structures, a strong correspondence with the features of the vDOS, of 0{u), and Pr{uj) 
of the bulk system can be found. Thus the chain and star are indeed representative of the 
low valence network, while the ring is characteristic of the network rich in three-functional 
junctions. 

As a consequence of the observed similarity between the chain and the X2 = 0.9 system we 
have the possibility to identify the different modes appearing in the vDOS. Figure |6] shows 
graphically the direction of the translational displacement of all particles belonging to the 
chain for characteristic frequencies that are located in different regions of the vDOS. One 
can identify the progression typical of one dimensional chains on increasing frequency, even 
if the purity of the modes is significantly altered by the curvilinear structure. Indeed, the 
graphical inspection of the displacements, as well as the specific values of Pr{u!) and 0{u)), 
consistently suggest a strong resemblance with a long wavelength acoustic transverse phonon 
(a), a short wavelength acoustic transverse phonon (b), a longitudinal acoustic phonon (c) 
and an optical-like longitudinal phonon (d). Of course, the curvilinear nature of the chain 
introduces some disorder that alter the pureness of the modes. 

Finally, we show in Fig. [7] three different modes of the ring, located in different frequency 
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FIG. 7. Three characteristic vibrational modes of the ring. 



regions. Apart from the low frequency mode which has a large participation ratio, the other 
modes describe deformations of the structure. The different values of 0{u) indicate if the 
deformation takes place along the plane of the ring or if it is mostly perpendicular to it. 

In order to further understand the origin of the peaks in the vDOS of our gel, it is useful 
to compare it to the one of a one-dimensional linear chain in the limit of very large A^. 
For this we have calculated the vDOS for a system of 1000 bi-functional particles located 
on a straight linear chain oriented in the x— direction. Neighbouring particles i and j are 
located at the distance corresponding to the minimum of V{i,j), where V{i,j) is the same 
inter-particle potential used in the simulation of the gel, see Eq. ([T]). The presence of angular 
degrees of freedom, which strongly couple to translational degrees of freedom, significantly 
alters the behavior of the vDOS with respect to the well-known text-book case of the one 
dimensional chain of particles interacting only via harmonic springs [1] . To evaluate the role 
of the coupling, we first start by studying the corresponding system in which the coupling is 
missing. This is realized by setting for the diagonalization the translation-orientation blocks 
of the Hessian to zero [H]. Beside the translational and rotational vDOS, we have also inves- 
tigated the dispersion relation, by analyzing the periodicity of the different eigenvectors. For 
this we have calculated the space Fourier transform of each eigenvector and have associated 
the position of the peak in A;— space to the corresponding eigenvalue u. In the case of trans- 
lation, motion along the chain direction provides information on the longitudinal branch, 
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FIG. 8. (a) Dispersion relation for the linear chain, calculated by estimating the spacial periodicity 
along the x direction (longitudinal) and along the y ox z (transverse) direction. The dispersion 
relations obtained by evaluating the periodicity for the two angles coincides with the the y ox z 
results, confirming the mixed nature of these modes. The inset shows the corresponding dispersion 
relation curves for the decoupled system in which we diagonalize only the 3A^ x 3A^ block of the 
Hessian associated to translational degrees of freedom as well as the 2N x 2N block associated to 
rotational degrees of freedom, (b) vDOS of a linear and of a bended chain of = 1000 particles. 
The vDOS of the network with X2 = 0.9 is also shown for comparison. The inset shows the vDOS 
calculated for the decoupled system. Note that the figure shows only the region u < 1.2 and 
therefore the van Hove singularity associated with the longitudinal modes and which is located at 
u ~ 4.4 is not visible. 

while displacements in the two orthogonal directions provide information on the transverse 
branches. Similarly, the periodicity of the two angular displacements of the particles along 
the chain provide information on the rotational modes. The dispersion relations and vDOS 
for the decoupled system are shown in the insets of Fig. |8^ and b. The longitudinal and 
transverse sound velocities vi and ft, calculated from the slopes of the dispersion curves at 
small k (or equivalently from the square root of the second derivative of the potential with 
respect to the chain direction x for the longitudinal modes and with respect to y oi z for 
the transverse ones) are found to be ^ 2.24 and Vt — 0.51. The significantly smaller value 
of ^a^a'"''' as compared to ^^y^ gives rise to a Vt that is much smaller than vi and hence 
a transverse dispersion relation located at frequencies much smaller than the longitudinal 
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one. This strong difference is plausible since it is much easier to deform the chain in its 
orthogonal direction than in its parallel direction. The translational vDOS of the system, 
inset of Fig. [8]d, shows a fiat density of states, followed by a van Hove singularity around 
w ~ 1 (the boundary zone of the transverse band) and around u ^ 4.5 (the boundary of the 
longitudinal band, not shown in the figure). The rotational modes have an optical character 
in that they have, at A; = 0, a finite frequency (inset of FigjS^). With increasing wave- vector 
the modes become softer, before they degenerate at the zone boundary. 

Next we come back to the original model, i.e. we include the coupling between translations 
and rotations when the Hessian is diagonalized. We see that this coupling induces strong 
changes in the dispersion relations, see main panel of Fig. |8^. The coupling makes that there 
are now two transverse optical modes that at A; = have the frequency of the rotational 
modes of the uncoupled system. Similarly the coupling makes that the rotational degrees 
of freedom can also have an acoustic branch, with a frequency at the zone boundary that is 
given by the u of the uncoupled system. In contrast to the rotational and transverse modes, 
no significant change in the longitudinal excitations are seen. As a result of the coupling 
between the various modes, the vDOS at low u becomes populated by modes which arise 
from the strong coupling between the low frequency transverse translational modes and the 
optical rotational modes (Fig. |8]d), thus explaining the presence of high intensity at low 
frequencies. 

We have also calculated the vDOS of a long bended chain, and the results are also 
reported in Fig. [Sjo. We see that the features of the linear chain are retained in the vDOS 
of the bended chain, even if the van-Hove singularities are now smeared out by the absence 
of linearity. Interestingly enough, the main features of the linear chain are also retained in 
the bulk system, whose broad bands can qualitatively associated to the linear chain bands, 
providing a strong support for the interpretation of the disordered network modes discussed 
previously (see Fig. |3]). 

IV. CONCLUSIONS 

In this article we have studied in detail the vDOS of a disordered system composed by 
poly-disperse chains of connected bi-functional particles, joined together via three-functional 
junctions. The relative concentration of bi- and three- functional particles makes it possible 
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to control the average distance / between junctions and explore the effect of the network 
mesh on the vDOS. 

We have found a clear evolution of the vDOS with /, with a progressive structuring of g{u)). 
To provide a characterization of the modes in different frequency regions, we have calculated 
the participation ratio and the mode amplitude in the bond direction, which has allowed 
us to distinguish between localized and delocalized modes and transverse and longitudinal 
excitations. We have found support for the resulting classification via the analysis of specific 
structures of the networks, i.e. chains, stars of chains and rings, finding that most of the 
features of the bulk vDOS can be associated with specific modes of these simpler structures. 

We have also found very helpful the comparison with the vDOS of an infinite chain of 
particles interacting with the same potential. The comparison between the vDOS calculated 
in the presence and in the absence of rotational-translational coupling helps clarifying the 
importance of such a coupling and the nature of the bands observed in the bulk disordered 
network. In particular, we observed that the coupling between translation and rotation 
significantly increases the vDOS at low u, which becomes populated by transverse modes, 
with frequencies which are strongly decreased (as compared to the uncoupled system) by 
the possibility of rotating the particles in coherence with the translational displacements. 
Interestingly, this band is clearly retained also in the vDOS of the bulk system. 

This study suggests that properties of the disordered systems can partially be traced 
back to properties of the corresponding crystal, with an adequate smearing introduced by 
the absence of long range periodicity and by the finite length of the chains. It also calls 
attention on the important role of the coupling between translation and rotation in the low 
frequency region of the spectrum, i.e. in the u region where the boson peak is observed in 
glasses. For the present model it is this coupling that leads to a strong intensity of the vDOS 
at low frequencies. However, this mechanism to generate low-frequency modes is certainly 
not the only one, since, i.e., disorder in the coupling constants or in the masses will also give 
rise to soft modes as has been shown for other systems [T71 HH] . 
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